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Recent measurements of the centrality and rapidity dependence of single inclusive jet produc¬ 
tion in p-fPb collisions at the LHC have revealed large and non-trivial nuclear modification of the 
production cross section for this process. In this paper, we explore to what extent such nuclear 
modification can be understood by the framework of standard cold nuclear matter effects, in par¬ 
ticular initial-state cold nuclear matter energy loss. We demonstrate quantitatively that theoretical 
calculations which include medium-induced radiative corrections can describe rather reasonably the 
attenuation of the jet production yields in the large transverse momentum region in d-|-Au collisions 
at RHIC and p-|-Pb collisions at the LHC for central to semi-central collisions. We further show that 
the observed scaling behavior of the nuclear modification factor as a function of the total jet energy 
Pt cosh( 3 /) for various rapidity intervals has a natural explanation in the picture of cold nuclear 
matter energy loss. On the other hand, the observed enhancement in peripheral collisions is not 
described in this picture and could have a different origin. 


I. INTRODUCTION 

High energy jet production in p-l-A collisions has long been regarded as a valuable baseline for studies of jet quenching 
and for constraining the properties of the hot dense medium created in A-l-A collisions at RHIC and LHC [j|, , since 

final-state effects associated with the quark-gluon plasma (QGP) are expected to be absent/suppressed. However, 
recent esroerimental results on the centrality dependence of high energy jet cross sections in p-|-Pb collisions at 
LHC d, 3] and in d-|-Au collisions at RHIC d, seem to challenge the expectation of negligible or small cross 
section modification. Highly nontrivial and large nuclear effects from different centrality selections are observed at all 
transverse momenta pt at forward (in the direction of the proton beam) rapidities and for large pt at mid-rapidity, 
and they are manifest as suppression of the jet yield in central events and enhancement in peripheral collisions Q • 

So far there is no theoretical approach which can successfully describe all aspects of these surprising results. In 
particular, they are challenging to explain within the framework of nuclear parton distribution functions (nPDFs) Q, 
which follows the usual leading-twist perturbative QCD factorization Q and attributes all nuclear effects to universal 
nPDFs in the large nucleus. There are also other attempts to understand the experimental data. For example, it was 
hypothesized 0 that the observed effects might arise from a suppression of the soft particle multiplicity in events 
with high energy jets, though the detailed mechanism of such suppression needs further detailed investigation. The 
correlation between the hard process rates and the soft event activity has also been studied in 0- Independently, 
it has also been argued [llljl3l| that parton configurations in the projectile proton containing a parton with large 
momentum fraction x interact with a nuclear target with a reduced cross section, which may give rise to the observed 
effects. At the same time, studies of kinematic bias on the centrality selection of jet events in p-l-A collisions at the 
LHC have been performed in [l3 |. Although the model is not based on any dynamical mechanism and fails in the 
case of peripheral collisions, it does capture the trend in the data at the LHC from central to semi-central collisions. 

In this paper, our goal is to explore to what extent the nuclear modification of jet production in p-l-A collisions 
observed at both RHIC and LHC can be described by the standard cold nuclear matter (CNM) effects. In p-l-A 
collisions, the production yield of single inclusive jets is affected by a number of CNM effects, all of which have clear 
physical/dynamical origin, mostly centered around the idea of multiple parton scattering. They include dynamical 
nuclear shadowing [T^ [T^ , Cronin effect [l3, [3 ’ incoherent multiple scattering mill and CNM energy loss [2ll- [ 2 ^ , 
which can all be calculated from the elastic, inelastic and coherent scattering processes of partons in large nuclei 30- 
l33| . It was found in our previous studies that dynamical shadowing, incoherent multiple scattering as well as Cronin 
effects play important roles only in the small and intermediate transverse momentum regions (32| . At very large 
transverse momentum pp, only CNM energy loss effects that follow the energy dependence obtained in [2^, [3 are 
relevant in modifying particle production yields in p-l-A collisions. Therefore, results on high pp jets in p-l-A collisions 
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should provide further important insights into the nature of CNM energy loss and differentiate between theoretical 
models. 

To this end, we study the phenomenological consequence of CNM energy loss effects on the high pT jet production 
in p+A reactions. We compare our theoretical calculations to sin^ inclusive jet production in d+Au collisions at 
RHIC at 200 GeV [l^ and in p+Pb collisions at LHC at 5.02 TeV Q from central to forward rapidity. We find that 
CNM energy loss leads to significant suppression of jet production in p+A collisions which is qualitatively consistent 
with both RHIC and LHC experimental data. In particular, the observed scaling behavior of the nuclear modification 
factor as a function of the total jet energy for various rapidity intervals finds its origin in the CNM energy loss picture. 
However, the observed enhancement in peripheral collisions is not described in this picture and needs to be further 
explored. 

The rest of our paper is organized as follows. In Sec. H, we first review the perturbative QCD formalism for single 
inclusive jet production in p+p collisions. We then discuss how the cold nuclear matter energy loss is implemented in 
this formalism. In Sec. HI, we present our phenomenological results of high pT jet production in d+Au collisions at 
RHIC and p+Pb collisions at the LHC. We demonstrate quantitatively to what extent the cold nuclear matter energy 
loss effects can explain the experimental data. We summarize our paper in Sec. IV. 


II. COLD NUCLEAR MATTER ENERGY LOSS IN P+A COLLISIONS 


To leading order (LO) in the framework of factorized perturbative QCD, single inclusive jet production in p+p 
collisions, p{Pa) + p{Pb) —t jet{y,pT) + X, can be written as [s^. 


dyfPpT 


= [ —fb/pixb,p)Hab^cis,i,u)S{s + t + u), 

S^JXa' JXb' 


( 1 ) 


where y and px are the rapidity and transverse momentum of the jet, respectively, j, represents the sum over all 
parton flavors, s = {Pa + Pb)'^ is the center-of-mass energy squared and s, t, u are standard Mandelstam variables 
defined at the partonic level. fa/p{x,fj,) are the parton distribution functions with p being the factorization scale. 
Hab^c are the short-distance hard-part functions for two partons of flavor a and b to reproduce a jet (at LO it is just 
the leading parton c), which are listed in Ref. [s^. In principle, a LO calculation of single inclusive jet production 
has certain limitations since only at next-to-leading order (NLO) the QCD structure of the jet starts to play a role 
in the theoretical description of physical observables , for example a LO calculation cannot be used to study 

the dependence of CNM effects on the jet cone radius R. Unlike final-state quark-gluon plasma effects, however, the 
dependence of initial-state effects on the jet cone radius is expected to be very small [^. Furthermore, quantitatively 
the NLO and LO cross sections are very similar for intermediate jet radii R ^ 0.4. For this reason LO calculations 
are very useful to discuss qualitatively the cold nuclear matter effects [i^. At the same time in the large pr region, 
it has been shown that the dependence of nuclear modihcation on the jet cone size is significantly reduced compared 
to that in the small and moderate px region. In this paper, we focus on the study of jets with very large transverse 
momentum px, and we thus take the advantage of the LO calculations for simplicity. 

In inclusive jet production in p+A collisions, the energetic partons before and after the hard collisions undergo 
multiple scattering and lose energy due to medium induced gluon bremsstrahlung. These processes are referred to 
as initial-state and final-state parton energy loss, respectively. A derivation to all orders in the correlations between 
the multiple scattering in nuclear matter of initial-state and final-state energy loss was performed in Ref. psl] using 
the Gyulassy-Levai-Vitev (GLV) reaction operator approach [4l|, see Fig. |T] for illustration of those cases. The 
bremsstrahlung pattern of Bertsch-Gunion beam jets, partons that only undergo soft interactions in nuclear matter, 
was also obtained. It was demonstrated that at high energies the fractional energy loss AE/E of beam jets is the 
largest. The resulting gluon bremsstrahlung may be partly or fully responsible for the soft particle production in p+A 
and A+A reactions and the observed azimuthal asymmetries [^, 1^ , however it is not relevant to high px particle 
and jet production. For large processes shown in Fig.[T] initial-state energy loss is much larger and dominates over 
the hnal-state energy loss when A —> oo. Note that in phenomenological applications of cold nnclear matter energy 
loss one must consider the large target nucleus rest frame, where the energies of the incident and outgoing parton are 
very high [2^ and, thus, final-state cold nuclear matter energy loss can indeed be neglected. Because of this, we will 
focus on the initial-state energy loss as it is central to our paper. 

The starting point of initial-state cold nuclear matter energy loss computation is the integral form of the double 
differential medium-induced gluon bremsstrahlung spectrum |25l | to first order in opacity: 


k+ 


d^Ns 
dk+Cpk I 


CrUs 


d^qx 


__ gj - gj_ • kx k+ 

l^g k^ikx - qx)"^ "k^ikl - qj_) k^Xc 


?eff 


gi 


-2 


■sm 


klL 

k+ 


( 2 ) 











3 





FIG. 1. Three distinct cases of medium-induced bremsstrahlung discussed in Ref. [2^ are illustrated. Top left: initial-state 
energy loss in the nucleus followed by a large process, resulting in the production of high-pr or high-S t particles and jets. 
Top right: final-state energy loss in the nucleus, after a hard collision. Bottom: the case of asymptotic on-shell Bertsh-Gunion 
beam jets that do not undergo a hard scattering process is also shown. 


where is the large light cone momentum of the radiated gluon, k±_ is the transverse momentum relative to the 
direction of the energetic parent parton, and g_L is the transverse momentum transfer between the propagating 
parton and the nuclear medium of size L and uniform density. Cp is the quadratic Casimir in the fundamental and 
adjoint representations of SU(3) for quarks and gluons, respectively. Xg is the gluon scattering length of 0(1 fm) 
and ^ is the typical transverse momentum transfer per interaction representing the interaction strength between the 
propagating parton and medium. On the other hand, the finite range of integration = 'C^jet/2 leads to 

+ Qmax)/Qmax [Hi- Effective parton mass due to interactions in matter is included as —>■ k‘j_+ 

The gluon number/intensity spectra and average gluon number/parton energy loss are then obtained as {k^ ~ 2a;): 
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In heavy ion collisions, where final-state energy loss effects dominate, phenomenological applications have improved 
upon the use of the mean energy loss AE in Eq. O- This approach can be generalized to initial-state energy loss. 
With the differential distribution of radiated gluons dN^/duj in Eq. ([3]) at hand, one can calculate the probability 
density Pq,g{d) for quarks and gluons to lose a fraction of their energy e = AEi/E due to multiple gluon emission 
in the Poisson approximation. We can further evaluate the mean energy loss fraction as follows: 

{eq,g) = imtiai-state ^ ^ ^ deePq,g{e), where ^ dePq,g{€) = 1 . (5) 

Note that in Eq. ([5]) with the subscripts q, g we explicitly indicate that quarks and gluons radiate different numbers 
of gluons and lose different fractions of their energy. This was implicit in Eqs. ([3]), (|3]). If the incident parton a (either 
quark or gluon) loses a fractional energy e, to satisfy the same final-state kinematics, it must carry a larger fraction 
of the colliding hadron momentum and, in turn, a larger value of momentum fraction Xa- This can be implemented 
in the cross section calculations in Eq. o in the following way [s^: 

fq/p{Xa,fJ.)^J^ dePq{e) , fg/p{Xa,g.) ^ d€Pg{€)fg/p(^-^^,H 


( 6 ) 
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Such an implementation has taken into account the fluctuation in the cold nuclear matter energy loss via Pq^g{e), and 
can be computationally demanding in the calculations of the nuclear modification factors in the next section. The 
result of energy loss fluctuations is a reduction in the attenuation of the particle and jet cross sections, effectively 
corresponding to a smaller fractional energy loss than that in Eq. To speed up the evaluation, we implement this 
effect as a momentum fraction shift in the PDFs, 


fq/p{Xa, fJ') fq/p 



fg/p{Xa, fg/p 



(7) 


where we follow Refs. [1^ use eq,g,eS = 0.7 {eq,g) with the mean energy loss fraction (eg,g) given by Eq. ([5]). 

This result immediately implies that the nuclear modification of single inclusive jet production in p+A collisions 
depends not only on the magnitude of initial-state cold nuclear matter energy loss, but also on the steepness of the 
parton distribution functions. In particular, large suppression is expected for jet production at forward rapidity and 
large transverse momentum, where the parton momentum fraction Xa in the proton is large and fq^g/p{xa, fJ-) is steeply 
falling. We will see these effects explicitly in the next section, where we evaluate the dependence of cold nuclear matter 
energy loss effects on the jet rapidity y and transverse momentum pT- 


III. NUMERICAL RESULTS AT RHIC AND LHC 


In this section we present numerical calculations of initial-state parton energy loss effects on inclusive jet production 
in p-|-A collisions and compare to experimental data from d-|-Au collisions at RHIC and p-fPb collisions at the LHC. 
To start, we define the usual nuclear modification factor RpA- 


1 daf^/dycPpT 

“ (iVcon) daf^Jdyd^pT ’ 


( 8 ) 



FIG. 2. Rep for inclusive jet production in d-|-Au collisions at ^/s = 200 GeV in central (left), mid-central (middle) and mid¬ 
peripheral (right) events. Data are from PHENIX [^. Yellow bands correspond to the theoretical uncertainty of initial-state 
parton energy loss in the range 0.175 < ^ < 0.7 GeV. 



FIG. 3. Same as in Fig. [21 but for i?dAu- 
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FIG. 4. Rep for inclusive jet production in p+Pb collisions at ^/s = 5.02 TeV in central (left), mid-central (middle) and 
mid-peripheral (right) events. Four different rapidity intervals (—0.3 < y < 0.3, 0.3 < y < 0.8, 0.8 < y < 1.2 and 1.2 < y < 2.1) 
are presented. Data are from ATLAS collaboration at the LHC [^. 


where (A^coii) is the average number of binary nucleon-nucleon collisions and is computed at a given centrality, and 
/dy(PpT and da^^^/dyd^pr represent the differential cross section of jet production in p-l-A and p+p collisions, 
respectively. The deviation of i?pA from unity reveals non-trivial nuclear effects in p-l-A collisions. Following experi¬ 
mental measurements, cold nuclear matter effects can also be quantified by the central-to-peripheral ratio Rep'- 


1 da?^/dyd'^pT 
Rcoii da?^^/dyd'^pT 


(9) 


where i?coU represents the ratio of A^^oii in a given centrality interval to that in the most peripheral interval. It is 
given by the expression Rcoii = (^con‘)/(^mn)- 

In both p-hp and p-l-A collisions, we use CTEQ6L1 parton distribution functions [1^. We choose the factorization 
and renormalization scales to be equal in the rest of this section and set them as p = px- In the numerical evaluation 
of initial-state parton energy loss, we use Ag = 1 fm for the gluon mean free path and vary the interaction strength 
between the propagating jet and the QCD medium by changing the typical momentum transfers ^ in the range of 
0.175 to 0.7 GeV. Such a range has been shown to give a good description for the nuclear modification of Drell-Yan 
dilepton production in fixed target experiments [^ . 

































































6 


O ATLAS 2.1 < y < 2.8 
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FIG. 5. Same as Fig. [4] but for the three other rapidity intervals (2.1 <y < 2.8, 2.8 < y < 3.6 and 3.6 <y < 4.4). 


d-l-Au 

Centrality 

0 - 20% 

20 - 40% 

40 - 60% 

60 - 88% 

(Neen) 

15.1 

10.2 

6.6 

3.2 

L [fm] 

10.2 

6.8 

4.1 

1.5 

p-fPb 

Centrality 

0 - 10% 

20 - 30% 

40 - 60% 

60 - 90% 

(Vcoll) 

14.6 

10.4 

6.4 

3.0 

L [fm] 

10.6 

7.3 

4.2 

1.6 


TABLE I. Estimated effective path length L for different centralities in d+An (top table) collisions at RHIC ^/s = 200 GeV 
and p+Pb (bottom table) collisions at the LHG ^/s = 5.02 TeV. (Vcoii) for RHIC and LHC are consistent with those from 
PHENIX collaboration at and ATLAS collaboration at 0], respectively. 


To calculate the centrality dependence of cold nuclear matter energy loss in p+A collisions, we use a Glauber 
model [U to determine the effective path length L in Eq. ([2]) for jet production at a given centrality class, in 
which we choose the inelastic nucleon-nucleon scattering cross section am = 42 mb at RHIC energy of y/s = 200 GeV 
and ain = 70 mb at the LHC energy of y/s = 5.02 TeV. The calculated effective path length L for the corresponding 
average number of binary collisions (iVcoii) are summarized in Table. HI for both d-|-Au collisions at RHIC and p-|-Pb 
collisions at the LHC. (Vcoii) for RHIC and LHC are taken from the PHENIX collaboration and the ATLAS 
collaboration Q, respectively. 

We are now ready to compare our theoretical results with experimental data. In Eig. we plot the central-to- 
peripheral ratio Rep for jet production in d-|-Au collisions at y/s = 200 GeV as a function of the jet transverse 
momentum pT, where the peripheral class is taken to be 60 — 88% centrality. The experimental data are from the 
PHENIX collaboration Q at central rapidity —0.3 < y < 0.3 in three different centrality intervals: 0 — 20% (left), 
20 — 40% (middle), 40 — 60% (right), respectively [i^. The yellow bands are our theoretical calculations, representing 
the uncertainty of typical momentum transfers between the propagating parton and nuclear medium in the range 
0.175 < ^ < 0.7 GeV as mentioned above. This momentum range has to be further constrained by unambiguous 
experimental measurements, e.g., the anticipated E906 Drell-Yan results from Eermilab fixed target experiments. 
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FIG. 6. i?ppb for inclusive jet production in p+Pb collisions at ^/s = 5.02 TeV in central (left) and minimum bias (right) 
events. Five different rapidity intervals (—0.3 < y < 0.3, 0.3 < y < 0.8, 0.8 < y < 1.2, 1.2 < y < 2.1 and 2.1 < y < 2.8) are 
presented. Data is from ATLAS collaboration at the LHC [j. 


As we have argued, cold nuclear matter energy loss leads to suppression of the jet production cross section. This 
attenuation becomes stronger in the most central collisions due to the increase in the effective path length L, as 
shown in Eq. ©■ Our numerical calculations give a reasonable description of the PHENIX data for i?cp. However, as 
shown in Fig. [31 it fails to describe the nuclear modification factor i?dAu for non-central collisions. In particular the 
enhancement in the peripheral collisions at 40 — 80% centrality class is not described in this picture. 

With the same set of gluon scattering length Xg and in-medium momentum transfers ^ in cold nuclear matter, 
we summarize our theoretical calculations of Rep in p-l-Pb collisions at LHC y/s = 5.02 TeV as a function of pr in 
Figs. [Handini where the peripheral class is taken to be 60 — 90% centrality. The experimental data are from ATLAS 
collaboration at the LHC 0. Seven rapidity intervals, —0.3 < y < 0.3, 0.3 < y < 0.8, 0.8 < y < 1.2, 1.2 < y < 2.1, 












































2.1 < y < 2.8, 2.8 < y < 3.6 and 3.6 < y < 4.4 in three different centrality bins (0 — 10%, 20 — 30%, 40 — 60%) 
are considered. Our theoretical calculations give a reasonable description for most of the data, and provide support 
for the cold nuclear matter energy loss picture. In particular, ATLAS data show a clear dependence of nuclear 
modification on the centrality, which seem to follow our expectation based on effective path length dependence in 
the cold nuclear matter energy loss formalism (^ . . In all centrality intervals the central-to-peripheral ratio i?cp 

decreases with increasing transverse momentum pr and rapidity y. Such dependences also find natural explanation 
in the cold nuclear matter energy loss picture. As we have emphasized in the last section, energy loss is very sensitive 
to the steepness of PDFs of the projectile proton. As pt and y increase, the sampled parton momentum fraction 
in the projectile proton Xa becomes larger, so does the the slope of the PDFs \dfa/pixa, p)/dxa\, and, thus, a larger 
suppression is expected. We notice here that in the highest px bins of the most central collisions in the most forward 
rapidity region, the nuclear suppression R^p can be as low as 0.2. This implies that initial-state cold nuclear matter 
effects in inclusive jet production may also be significant in A-|-A reactions in this kinematic region (40l| . 

We further plot the nuclear modification factor i?ppb in the most central (0 — 10%) (left) as well as minimum 
bias (right) collisions in Fig. [6l The same patterns as in i?cp and reasonable agreements between our theoretical 
calculations and experimental data are found. The agreement is best in central collisions, shown in the left panels 
of Fig. |6l A nuclear PDF-based calculation does not reproduce this suppression. However, similar to RHIC findings, 
anomalous nuclear enhancement of single jet production in peripheral collisions is also observed at the LHC. Such an 
enhancement is again beyond what we can expect from either the initial-state parton energy loss picture or the nPDF 
model. In minimum bias collisions, shown in the right panels of Fig. [51 there is no clear indication that i?ppb deviates 
form unity, as expected in energy loss calculations. This is consistent with the nPDF large behavior. Within 
statistical and systematic uncertainties, the data is, however, also compatible with the upper edge of the CNM energy 
loss band. 
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O ATLAS 0.3 < y < 0.8 
ATLAS -0.3 < y < 0.3 


500 1000 1500 2000 200 400 600 800 1000 

Pj*cosh(y) (GeV) 


QL 


ATLAS 2.1 < y < 2.8 
o ATLAS 1.2 < y < 2.1 
ATLAS 0.8 < y < 1.2 
o ATLAS 0.3 < y < 0.8 



000 1500 2000 

PTAcosh(y) (GeV) 


FIG. 7. Rep (top) and Rppb (bottom) in 0 — 10% p-l-Pb collisions as a function of pr x cosh(y) are shown at different rapidities 
in the most central collisions (0 — 10%). Seven rapidity intervals are presented for Rep and four for Appb. Data is from ATLAS 
collaboration at the LHC [j. 
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In search of a universal relationship between the nuclear modification factors in different rapidity ranges ATLAS 
has further plotted Rep and i?pPb as a function of the jet kinematic variable pt x cosh(j/) [^, which is approximately 
the total energy of the jet. They find that, when plotted against this variable, the nuclear modification factors in the 
central and forward rapidity regions fall along the same curve. However, Rep and i?ppb at backward rapidity do not 
follow this trend. These patterns suggest that the observed nuclear modification may depend on the initial parton 
kinematics, especially the parton momentum fraction Xa in the projectile proton. Such a behavior is exactly what 
one can expect in the initial-state parton energy loss picture, which we will now elaborate in more detail. From the 
partonic kinematics at LO one can easily derive in the forward rapidity region, where Xa ^ Xb, 

Pt X cosh( 2 /) oc Xa\/s , (10) 

and we have used the fact that the cross section is dominated by configurations where the observed and recoil jets are 
at approximately equal rapidities. Therefore, the kinematic variable pt x cosh(y) introduced by ATLAS is simply the 
projectile parton momentum fraction Xa times a constant factor. Since initial-state parton energy loss is equivalent 
to a shift in Xa, as shown in Eq. 0, the nuclear modification factors will follow the same trend once they are plotted 
as a function of Xa- This is precisely what we have observed in Fig. [71 where we show the comparison between our 
theoretical calculations of Rep and i?ppb and the experimental data, against the quantity px x cosh(j/). As one can 
see, all theoretical calculations for different rapidity intervals are on top of each other, following the correct scaling 
with Xa- There are clearly some quantitative difference between our theoretical calculations and experimental data, 
especially in Rep, which we attribute partly to the unexplained anomalous enhancement of jet production in peripheral 
collisions. With this caveat, the reasonable agreement on both the trend and the size of the nuclear modification in 
p-|-A collisions is still very encouraging. On the other hand, in the backward rapidity region Xb is comparable to (even 
larger than) Xa, and one does not have a simple proportionality between Xa and the total jet energy. Thus, a scaling 
behavior, as observed in forward rapidity region, should not be expected. This is exactly what ATLAS collaboration 
has observed. 


IV. SUMMARY 


In this paper we studied to what extent the nuclear modification of inclusive jet production in p-l-A collisions can be 
described by standard cold nuclear matter effects. In particular, for jet production at large transverse momentum px 
we evaluated the nuclear modification factor i?pA and the central-to-peripheral ratios Rep of inclusive jet production 
due to initial-state cold nuclear matter energy loss in d-|-Au collisions at RHIC and p-fPb collisions at the LHC. We 
found that our theoretical calculations can capture quantitatively the bulk of the observed modifications for both 
RHIC and LHC experiments from central to semi-central collisions. The agreement is good in central reactions and 
the upper edge of the CNM energy loss band is still consistent with data in minimum bias reactions if the statistical 
and systematic uncertainties are taken into account. On the other hand, a nPDF-based calculation is consistent with 
the minimum bias data in inclusive jet production form RHIC and LHC, but fails to describe central collisions. What 
is particularly important to note is that the observed scaling patterns of Rep and i?pPb against the total energy of 
the jet px X cosh(?/) oc Xa (the parton momentum fraction in the projectile proton) in the forward rapidity region in 
p-|-Pb collisions find their natural explanation in the picture of cold nuclear matter energy loss. We further note that 
Xa ~ xp, Feynman x, at forward rapidity and similar scaling of the nuclear attenuation is expected for large xp- It 
has indeed been observed for different final states jsoj i . The encouraging comparison between theoretical calculations 
and experimental data indicates the significance of cold nuclear matter energy loss for understanding particle and jet 
production in p-|-A collisions. 

On the other hand, the observed nuclear enhancement in peripheral collisions is difficult to understand in our 
picture or the nPDF model. Such an enhancement might have a different origin, for example from “centrality bias”, 
and needs to be explored further. So far it has no first-principles explanation at the quantitative level. The anomalous 
enhancement of jet production also affects the central-to-peripheral ratio and is partly responsible for the small values 
of Rep- Experimentally, future results from RHIC, LHC and the fixed target p-l-A program at Fermilab will be 
essential to further elucidate the non-trivial nuclear modification of high pr/high invariant mass particles and jets in 
p-|-A reactions and to constrain the transport properties of cold nuclear matter in the current energy loss formalisms. 
Theoretically, it will be important to understand if there is centrality selection bias and what is its dynamical origin. 
As far as inelastic initial-state interactions are concerned, a step forward will be to go beyond the soft-gluon energy 
loss approximation and obtain the full medium-induced splitting kernels [5ll| . With such cold nuclear matter splitting 
kernels at hand, one can treat vacuum and in-medium parton showers on the same footing, following the progress 
recently made on the implementation of final-state QGP effects We leave such developments for future work. 
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